June, 2012 



SU(3) symmetry breaking and CP violation in D — > PP decays 

Hai-Yang Cheng 

Institute of Physics, Academia Sinica, Taipei, Taiwan 11529, ROC 

Cheng- Wei Chiang 

Department of Physics and Center for Mathematics and Theoretical Physics, 
National Central University, Chungli, Taiwan 32001, ROC 
Institute of Physics, Academia Sinica, Taipei, Taiwan 11529, ROC and 
Physics Division, National Center for Theoretical Sciences, Hsinchu, Taiwan 30013, ROC 

Abstract 

Evidence of CP violation in the charm sector has been observed recently by the LHCb and CDF 
Collaborations. Adopting the topological diagram approach, we study flavor SU(3) symmetry 
breaking effects in the weak decay tree amplitudes of singly Cabibbo-suppressed D — > PP decays. 
The symmetry breaking in the color-allowed and color-suppressed amplitudes is estimated with the 
help of the factorization ansatz, while that in the VF-exchange amplitude is done by fitting to related 
branching fraction data. We find that the W-exchange amplitudes stay in the second quadrant 
relative to the color-allowed tree amplitude, albeit there are two possibilities for one type of W- 
exchange amplitude. The weak decay penguin amplitudes, on the other hand, are evaluated within 
the framework of QCD factorization. Using the input of topological tree amplitudes extracted from 
the Cabibbo-favored decay modes and the perturbative results for QCD penguin amplitudes, we 
make predictions for the branching fractions and CP asymmetries of singly Cabibbo-suppressed 
modes. The predictions of branching fractions are generally improved from those in the SU(3) 
limit. We conclude that the direct CP asymmetry difference between D° — > K + K~ and D° — > 
7T+7T- is about -(0.139 ± 0.004)% and -(0.151 ± 0.004)% for the two solutions of lU-exchange 
amplitudes, respectively. We also find that the CP asymmetry of D° — > K°K° dominated by 
the interference between lU-exchange amplitudes ranges from -0.62 x 10" 3 to -1.82 x 10" 3 . We 
study phenomenological implications of two new physics scenarios for explaining the observed 
CP asymmetry in the charm sector, one with large penguin amplitudes and the other with a large 
chromomagnetic dipole operator. We find that the two scenarios can be discriminated by the 
measurements of CP asymmetries of a set of decay modes. 
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I. INTRODUCTION 



Recently, the LHCb Collaboration has reported a result of a nonzero value for the difference 
between the time-integrated CP asymmetries of the decays D° — > K + K~ and D° — > 7t + tt~ [l| 

AA C P = Acp(K + K-)-A C p{7r + 7T~) = -(0.82 ±0.21 ±0.11)% (LHCb) (1) 

based on 0.62 fb _1 of 2011 data. The significance of the measured deviation from zero is 3.5o\ This 
first evidence of CP violation in the charm sector was later confirmed by the CDF Collaboration 
based on a data sample corresponding to the integrated luminosity of 9.7 fb^ 1 Q| 

AA CP = -(0.62 ±0.21 ±0.10)% (CDF). (2) 
The time-integrated asymmetry can be written to first order as 

A C p(f) = af P (f) + { ^Sp , (3) 

where a^p is the direct CP asymmetry, a™p is the indirect CP asymmetry, (t) is the average decay 
time in the sample, and r is the D° lifetime. The combination of the LHCb, CDF, BaBar and 
Belle measurements yields a£p = -(0.025 ± 0.231)% and Aa$ P = -(0.656 ± 0.154)% Q|. 

It is important to explore whether the first evidence of CP violation in the charm sector is 
consistent with the standard model (SM) or implies new physics (NP). For some early and recent 
theoretical investigations, see Refs. [4426 ] . 

It is known that a reliable theoretical description of the underlying mechanism for exclusive 
hadronic D decays based on quantum chromo dynamics (QCD) is still not yet available. This is 
because the mass of the charm quark, being of order 1.5 GeV, is not heavy enough to allow for a 
sensible heavy quark expansion. Indeed, it does not make too much sense to generalize the QCD 
factorization (QCDF) and perturbative QCD (pQCD) approaches to charm decays as the l/m c 
power corrections are so large that the heavy quark expansion is beyond control. 

It turns out that a more suitable framework for the analysis of hadronic charmed meson decays 
is the so-called topological diagram approach, a powerful tool for a model-independent analysis. In 
this approach, the topological diagrams are classified according to the topologies in the flavor flow of 
weak decay diagrams, with all strong interaction effects included implicitly. Based on flavor SU(3) 
symmetry, this model-independent analysis enables us to extract the topological amplitudes by 
fitting to available data, to probe the relative importance of different underlying decay mechanisms, 
and to relate one process to another at the topological amplitude level. 

The salient point of the topological diagram approach is that the magnitude and the relative 
strong phase of each individual topological tree amplitude in charm decays can be extracted from 
the data. This allows us to calculate CP asymmetry at tree level in a reliable way, granting us an 
idea about the size of CP violation in charmed meson decays. Based on this approach, we have 
studied AAcp in D — > PP and D — > VP decays and obtained Aa^p « —0.14% and an upper 
bound around —0.25% in the SM [3]. A very similar result of Aa^p = —0.118% based on a 
variant of the diagrammatic approach was obtained recently in 22]. 

While many NP models 0-Q Jl3fl, H 

had been proposed to explain the measurement 



of large AAqp, it was argued in [1CJ, [l4|, [l6|, |2l|, |23j, |24| that large CP asymmetries in singly 
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Cabibbo-suppressed (SCS) D decays were allowed in the SM and the measured Aa^p could be 
accommodated or marginally achieved. In particular, it was advocated by Brod et al. 2J] the 
possibility of a large penguin amplitude in the SM. Denote d- and s-quark penguin contractions 
by Pd and P s , respectively. Under the assumption of large enhancement in P^ s relative to the tree 
amplitude, Brod et al. claimed that the sum of Pd and P s could explain Aa^p, while the difference 
between Pd and P s explains the large disparity in the rates of D° — > K + K~ and D° — > 7r + 7r _ . This 
would require that (Pd — Ps) be of the same order as the tree amplitude. Anyway, the interpretation 
of the seemingly large SU(3) breaking in the ratio R = T(D° — > K + K~)/T(D° — > ir + ir~) in terms 
of the difference of penguin contractions (Pd — P s ) is at odds with the common wisdom that 
the overall large SU(3) symmetry violation arises from the accumulation of several small SU(3) 
breaking effects in the T and E amplitudes. In this work we shall carefully examine the large 
penguin interpretation of R. 

As for the NP explanation of Aa^p, the authors in 18] argued that a large chromomagnetic 
dipole operator could be the best NP candidate to explain the data while satisfying most flavor 
physics constraints at the same time. In the present work we shall focus on the aforementioned 
pictures of large penguins and large chromomagnetic operator in SCS D — > PP decays. In par- 
ticular, we will investigate their phenomenological consequences in the CP asymmetries of these 
charmed meson decay modes, seeing if there are discernible differences in the two scenarios. 

Based on the topological diagram approach, we have recently studied direct CP asymmetries in 
the SM for SCS D — > PP and D — > VP decays [3| . In this work we will improve the previous anal- 
ysis by extracting the VF-exchange amplitudes directly from the decays D° — > K + K~ ,tt + it~ ,7t tt° 
and K°K° and considering explicitly all SU(3) breaking effects in SCS decay amplitudes. 

The layout of the present paper is as follows. After a brief review of the diagrammatic approach, 
we study various mechanisms responsible for the large SU(3) violation in the branching fraction 
ratio of D° — > K + K~ to D° — > ir + ir~ and fix the weak annihilation amplitudes in Section ITT1 
Unlike the topological tree amplitudes which can be extracted from the data, penguin amplitudes 
in this work are evaluated in the framework of QCD factorization as illustrated in Section IIIH We 
then discuss direct CP violation in SCS D — > PP decays in Section [IV] and possible explanations 
of the LHCb and CDF measurements in terms of new physics in Section |VJ Finally, Section IVII 
comes to our conclusions. 



II. DIAGRAMMATIC APPROACH 

A. Topological amplitudes 

It has been established sometime ago that a least model-dependent analysis of heavy meson 
decays can be carried out in the so-called quark diagram (or topological diagram) approach 
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291 ]. In this diagrammatic scenario, the topological diagrams can be classified into three distinct 



groups (see Fig. 1 of [15]). For our purposes, it suffices to consider two of them (see [3fjj for details): 
1. Tree and penguin amplitudes: color-allowed tree amplitude T; color-suppressed tree am- 
plitude C; QCD-penguin amplitude P; singlet QCD-penguin amplitude S involving flavor SU(3)- 
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singlet mesons; color- favored electroweak-penguin (EW-penguin) amplitude Pew', color-suppressed 
EW-penguin amplitude P^w- 

2. Weak annihilation amplitudes: VF-exchange amplitude E; W-annihilation amplitude A; 
QCD-penguin exchange amplitude PE; QCD-penguin annihilation amplitude PA; EW-penguin 
exchange amplitude PEew', EW-penguin annihilation amplitude PAew- 

It should be stressed that these diagrams are classified purely according to the topologies of 
weak interactions and flavor flows with all strong interaction effects encoded, and hence they are not 
Feynman graphs. In other words, all quark graphs used in this approach are topological and meant 
to include strong interactions to all orders, with gluon lines and quark loops included implicitly in 
all possible ways. Therefore, analyses of topological graphs can provide information on final-state 
interactions (FSI's). 

The topological amplitudes T, C, E, A are extracted from the Cabibbo-favored (CF) D — > PP 



decays to be (in units of 10 6 GeV) |3l(] (see also |32j]) 



T = 3.14 ± 0.06, C = (2.61 ± 0.08) e^ 152 * 1 )" , 

E = (1.53l°;0 8 7 ) e *( 122±2 )° , A = (0.39±°;^) e^ 31 -^ (4) 



for 6 = 40.4° 



331 ] . where cp is the r) — rf mixing angle defined in the flavor basis 



rf 




(5) 



with r] q = -^(uu + dd) and r] s = ss. The fitted \ 2 value is 0.29 per degree of freedom. 

The topological amplitudes C and E given in Eq. @ extracted from the data are much larger 
than those expected from naive factorization. In the factorization approach, the color-suppressed 
tree amplitude C has the expression 

C = ^a 2 (K7r) f K (m 2 D - ml)F^{ml). (6) 

Using the form factor for D to ir transition determined from measurements of semileptonic D meson 



decays to the tt meson 34J], we find 02(^71") = (0.82±0.02)e~*( 152±1 )° 31], to be compared with 02 = 
C2 + ci/3 ~ —0.11 in naive factorization. Likewise, weak annihilation diagrams should be helicity 
suppressed, whereas data imply larger sizes for them. This is because they receive large l/?n c power 
corrections from FSI's and large nonfactorizable contributions for For example, the topological 
amplitude E receives contributions from the tree amplitude T via final-state rescattering with 
nearby resonance effects. The large magnitude and phase of weak annihilation can be quantitatively 



and qualitatively understood as elaborated on in Refs. [371. l38| 



B. Flavor SU(3) symmetry breaking 

Under the flavor SU(3) symmetry, one can use the topological amplitudes extracted from the 
CF modes to predict the rates for the SCS and doubly Cabibbo-suppressed (DCS) decays. It is 
known that while the agreement with experiment for the branching fractions of SCS D — > PP 
decays is generally good (see the second column of Table Hill . there exist significant SU(3) breaking 
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in some of the SCS modes in the flavor SU(3) symmetry limit. For example, the predicted rates for 
the 7r + 7r _ and 7r°7r° modes are too large, while those for the K + K~ , ir + r]^ and K + t]^ modes are 
too small compared to data. The decay D° — > K°K° is almost prohibited in the SU(3) symmetry 
limit, but the measured rate is of the same order of magnitude as that of D° — > 7r°7r°. 

The most noticeable example of SU(3) breaking lies in the decays D° — > K + K~ and D — > 
7r + 7r~. Experimentally, the branching fraction of D° — > K + K~ is larger than that of D° — > it + it~ 



by a factor of 2.8 [35|, while their rates should be the same in the SU(3) limit. This is a long- 
standing puzzle since SU(3) symmetry is expected to be broken merely at the level of 30%. Without 
the inclusion of SU(3) breaking effects in the topological amplitudes, we see from Table ITU that the 
predicted rate of K + K~ is even smaller than that of tt + tt^ due to less phase space available to the 
former. 

The conventional wisdom for solving the above-mentioned long-standing puzzle is that the 
overall seemingly large SU(3) symmetry violation arises from the accumulation of several small 



and nominal SU(3) breaking effects in the tree amplitudes T and E (see e.g. |36l]). From the recent 
measurement of AAcp by LHCb and CDF, we learn that penguin diagrams in SCS decay channels 
do play a crucial role for CP violation. This leads some authors to conjecture that penguins may 
explain the rate disparity between D° — > K + K~ and D° — > 7r + 7r~ . 
To begin with, we write 

A{D° -> 7T+7T-) = \ d (T + E + P d + PE d + PA d ) nn + X S {P S + PE S + PA S )^ 

= ^(X d - A S )(T + E + AP)™ - ^X b (T + E + SP) W7r , (7) 

where X p = V* p V up {p = d, s,b) and 

AP = (P d + PE d + PA d ) - (P s + PE S + PA S ) , 

SP ee (P d + PE d + PA d ) + (P s + PE S + PA S ) . (8) 

Likewise, 

A(D° -> K+K~) = X d (P d + PE d + PA d ) KK + X S {T + E + P S + PE S + PA S ) KK 

= \{K~ X d )(T + E- AP) KK - l -X b (T + E + Y,P) KK . (9) 



The quantities AP and SP are denoted by Vj and 2Pj, respectively, in 2J]. Recently, an interest- 



ing picture has been proposed in [2j] that AP dominated by the difference of d- and s-quark penguin 
contractions of 4-quark tree operators can explain the large rate difference between D° — > tt + it~ 
and D° — > K + K~, while ^SP, the average of d- and s-quark penguin contractions, can account 
for the large direct CP violation Aa^p. Here we would like to examine the quantity AP. 

First, assuming SU(3) symmetry for the amplitudes T and E for the moment, a fit to the 
measured branching fractions B(D° -> vr+vr-) = (1.400 ± 0.026) x 10~ 3 and B(D° K + K~) = 



(3.96 ± 0.08) x 10" 3 [35( yields (in units of 10" 6 GeV) 

AP= 1.54e^ 202 °. (10) 

Therefore, |AP/T| ~ 0.5. Since the real part of AP is negative, it is clear that AP contributes 
constructively to K + K~ and destructively to 7r + 7r~. Next, assume that [/-symmetry breaking in 
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the amplitudes T + E follows the pattern 24 ] 



(T + E) W7T = (T + E)(l + le T ), (T + E) KK = (T + E)(l — \e T ), (11) 



where e T is a complex parameter with |e T | £ (0,0.3). It was shown in [24J that the relation 
|AP/T| ~ 0.5 still holds roughly. Hence, the large rate disparity between D° — > 7r + 7r~ and 
D° — > K + K~ is ascribed to the difference of d- and s-quark penguin contractions. In doing so, 
AP needs to be of the same order of magnitude as T. 

However, SU(3) breaking in reality does not necessarily follow the pattern exhibited in Eq. (jlip . 
In the factorization approach, various topological T amplitudes have the expressions 

T K7r = | fll (ft)AK-m|)^K), 

T KK = ^a x {KK)j K (ml-m\)F» K {m\\ (12) 

Ttttt = -j=ai(TTTr) f n (m 2 D -ml)F^ n (ml), 

for CF D° K--K+, SCS D° K + K~ and D° vr+vr- decays, respectively. Hence, 
T KK _ ai (KK) f K F DK (m 2 K ) m 2 D - m\ _ ^ a^KK) 

T T7T oiCtttt) U F^ijnl) m 2 D - ml ' a^mr) ' 1 ' 

where we have used the form- factor q 2 dependence determined experimentally from Ref. 34|. The 



parameter a\ has the general expression (see e.g. [481 ]) 



ai(PiP 2 ) = ( c i + §-) + §- ^[^(^) + ^-^(PM] + 0(l/m c ) corrections , (14) 
where V\ and H\ are vertex corrections and hard spectator interactions, respectively. In the 



diagrammatic approach, the parameter a\ can be extracted to be |ai(ET7r)| = 1.22 ± 0.02 31] 
from the data of CF D — > Ktt decays. Recall that c\ + y ~ 1.274 at the scale /i = m c ~ 1.3 
GeV. Therefore, it is evident that nonfactorizable contributions and l/m c corrections to a\ are 
rather small. It is thus reasonable to take the ratio a\{K K) / a^TTir) to be in the vicinity of unity. 
Neglecting SU(3) breaking in the VF-exchange amplitudes for the moment {i.e., E KK = E nn ), we 
get 1 

AP = 0.49e- il29 ° . (15) 

Therefore, |AP/T| = 0(0.15) rather than 0(1) in the presence of SU(3) violation in T amplitudes 
suggested by the factorization approach. 

Furthermore, SU(3) symmetry should be also broken in the VF-exchange and penguin annihi- 
lation amplitudes. This can be seen from the observation of the decay D° — > K !^ whose decay 
amplitude is given by 

A(D° -> K°K°) = X d (E d + 2PA d ) + \ S (E S + 2PA S ) , (16) 



1 In the same spirit, Bhattacharya ct al. Q found "P" = \ d (P + PA) d + X S (P + PA) S = (0.044 ± 0.023) + 
i(0.141 ± 0.036) after a fit to the rates of D° K + K~, tt+tt" and 7r°7r°. Then AP is obtained through 
the approximation AP sa "P"/A<2 « — "P"/ sin^c- 
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where E q refers to the VF-exchange amplitude associated with cu -4- qq (q = d,s). In the SU(3) 
limit, the decay amplitude is proportional to A;, and hence its rate is negligibly small, while ex- 



perimentally B(D° -4 K°K ) = (0.346 ± 0.058) x 10~ 3 |35fl. This implies sizable SU(3) symmetry 
violation in the VF-exchange and QCD-penguin annihilation amplitudes. Since the theoretical es- 
timate of AP is small [see Eq. (|40p below], we shall neglect AP and assume that the T and E 
amplitudes are responsible for the SU(3) symmetry breaking. Noting that E nn = Ed and E KK = E s 
in our notation, we can fix the VT-exchange amplitudes from the following four modes: K + K~ , 
7r + 7r~, 7r°7r° and K°K°. 2 Neglecting AP and A;, terms (see Table U for the topological amplitudes 
of the above four decay modes), a fit to the data yields two possible solutions 

(I) E d = 1.19e il5 -°°P, E s = 0.58e- il4 - 7 °P , 

(II) E d = lA9e il5 -°°E, E s = 1.62e- i9 - 8 °E . (17) 

The corresponding \ 2 vanishes as these two solutions can be obtained exactly. It has been noticed 
that a significant phase difference between Ed and E s is needed in order to fit the data of D° — > 



K !^ 39 1 and to account for the large rate difference between D° — > K + K and D° — > it + it 



31 ] . When comparing the predicted branching fractions of D° — > rjrj and D° — > r\r( (see Table [II]) , 



it appears that Solution I is slig htly more preferable, but Solution II is equally acceptable. 



We have argued in Ref. 3lf] that the long-distance resonant contribution through the nearby 
state /o(1710) could account for SU(3)-breaking effects in the ^-exchange topology. This has to 
do with the dominance of the scalar glueball content of /o(1710) and the chiral suppression effect 
in the ratio r(/ (1710) -4 7T7f)/r(/ (1710) -4 KK). 

To summarize, if SU(3) symmetry in T and E amplitudes holds or is broken in the manner 
depicted by Eq. (TIT]) , then a sizable AP of the same order of magnitude as T is needed to explain 
the data. If SU(3) violation due to decay constants, meson masses and form factors is taken into 
account in T amplitudes so that T KK /T n7T ~ 1.32 , it leads to |AP/T| ~ 0.15 . Finally, if the large 
rate difference between D° — > K + K~ and D° — > tt + -k~ is entirely accounted for by SU(3) violation 
in T and E amplitudes, AP will be negligibly small. Owing to the observation of D° — > K°K° 
through VF-exchange and penguin annihilation diagrams and the smallness of AP theoretically, we 
shall argue in this work that the last scenario is preferred. 

To estimate the effects of SU(3) symmetry violation in T and C amplitudes, we shall rely on 
the factorization ansatz. In this approach, the topological amplitudes T and C extracted from CF 
D — > Kn decays have the expressions 

T = ^(Kn) U(m 2 D - m 2 K )F DK (ml), (18) 



and Eq. ([6]) for the amplitude C. In [3j| we have illustrated SU(3) breaking effects in some selective 



SCS modes. For example, the relevant factorizable amplitudes for D + — > n + r]^ decays are 



2 However, the W-exchange contribution is missing in the topological amplitude expression of the D° 
K°K° decaj 
diagram PA 



K°K° decay in 16 1 as SU(3) symmetry is assumed for ^-exchange but not for the penguin annihilation 
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TABLE I: Topological amplitudes for singly Cabibbo-suppressed decays of charmed mesons to two 
pseudoscalar mesons where flavor SU(3) symmetry breaking effects are included. Summation over 
p = d, s is understood. For simplicity, flavor-singlet QCD penguin, flavor-singlet weak annihilation 
and electroweak penguin annihilation amplitudes have been neglected. 
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Cirri, = ^a 2 fg(rn 2 D -ml)F^ )7T (rn 2 r)q ), 

C^ a = ^a 2 f s (m 2 D -ml)F D *(ml), (19) 

where f q , f s are the decay constants of r\ q and rj s , respectively. We shall use the parameters f q = f n , 
f s = 1.325f n , m Vq = 741 MeV and ///„, = 783 MeV 40], and assume the form factor F^ 9 to be 



the same as F^ n . As pointed out in [3J, one needs SU(3) violation in weak annihilation to get a 
better agreement with experiment. For this purpose, we rely on the decay constants involved in 
the processes to estimate the SU(3) breaking effects in the ^-annihilation amplitude A. In CF 
D —7- PP decays, the topological amplitude A is extracted from Df — > K + K° which involves the 
decay constants f Ds and f K - For D + ->■ tt + t], we thus have A nri = (fD/fD s )(f-7rf q /fK) A = 115^4 
and, likewise, A^rf = ( fn/ fn„)( f-n-fx/ fx) A = 1.56A, where use of world averages fo = 213 MeV 



and Jd s = 248 MeV 4lf] has been made. Finally, the decay amplitudes read (see also Tabled]) 



A(D + ^vr+r/) = -=1^1^(0.816 T + 0.930 C + 1.15A) cos 0- V* S V US 1.285C sin 0, 
v2 

A{D + -> ir+rf) = -=V^ l V ud (0.8WT + 0.930 C + IMA) sin <p + V* s V US 1.285C cos </>. (20) 
v2 

From Table HH we see that the discrepancy between theory and experiment for B(D + — > Tr + r]^) is 



TABLE II: Branching fractions (in units of 10 -3 ) of singly Cabibbo-suppressed D — > PP decays. 
The column denoted by £> su(3) shows the predictions based on our best-fitted results in Eq. @ with 
exact flavor SU(3) symmetry, while SU(3) symmetry breaking effects are taken into account in the 
column denoted by £> su(3 )_ brcaking - The first (second) entry in D° — > r]r], rjrf, K + K~ and K°K° 
modes is for Solution I (II) of and E s in Eq. (|17p . Experimental results of branching fractions 
are taken from Ref. [351 ] . 
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4.62 ±0.31 


4.67 ±0.29 


D+ - 




5 


46 ±0.55 


8.93 ±0.85 


5.66 ±0.32 


Df- 




2 


85 ±0.28 


3.65 ±0.33 


2.42 ±0.16 


Df- 







73 ± 0.09 


0.73 ± 0.09 


0.62 ±0.21 


Df- 


-> K+T) 





79 ± 0.08 


0.84 ± 0.08 


1.75 ±0.35 


Df~ 


-> K+n' 


1 


02 ±0.17 


1.25 ±0.20 


1.8 ±0.6 



greatly improved. 

SU(3) breaking effects in the topological amplitudes for SCS D — > PP decays are summarized 
in Tabled] Theory predicted and measured branching fractions are given in Table HT1 3 While the 
agreement with experiment is improved for most of the SCS modes after taking into account SU(3) 
breaking effects in decay amplitudes, there are a few exceptions. For example, the predicted rate 
for D + — > K + K° becomes even worse compared to the prediction based on SU(3) symmetry. It 
is possible that the effective parameter a±(D + — > K + K°) is smaller than ai(D° — > Kir). Finally, 
we note in passing that the central values of the predicted branching fraction in the SU(3) limit 



given in Table [TI] are sometimes slightly different from those given in Table II of Ref. [3lJ . This is 



Our previous result of ai£ t) (D+ -> K+rj) is erroneous and it has been corrected in the erratum of 
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because the current work adds penguin contributions into the analysis. The errors associated with 
these branching fractions are larger than before because a Monte Carlo sampling method is used 
here for error estimation, instead of the simple error propagation method used in Ref. 3l[. The 
same Monte Carlo method is used throughout this paper for error estimation. 



III. PENGUIN AMPLITUDES IN QCD FACTORIZATION 

Although the topological tree amplitudes T, C, E and A for hadronic D decays can be extracted 
from the data, we still need information on penguin amplitudes (QCD penguin, penguin exchange 
and penguin annihilation) in order to estimate CP violation in the SCS decays. Unlike the tree 
amplitudes, it is more difficult to extract the topological penguin amplitudes reliably from the data. 
This is because the use of the topological approach relies heavily on SU(3) symmetry which leads 
to negligible penguin amplitudes in D decays. Consequently, the extraction of penguin amplitudes 
depends on SU(3) breaking effects in tree amplitudes. Indeed, we have shown in Sec. Ill Bl that the 
difference in penguin contractions AP is sensitive to how SU(3) symmetry breaking is treated in the 
tree amplitudes T and E. Therefore, we shall rely on theory to estimate the penguin contribution. 

With the advent of heavy quark effective theory, it is known that nonleptonic B decays can be 
analyzed systematically within the QCD framework. There are three popular approaches available 



in this regard: QCD factorization (QCDF) [42j|, perturbative QCD (pQCD) [43J and soft-collinear 



effective theory (SCET) 44f|. QCDF is phenomenologically quite successful in describing charmless 



hadronic B decays (see Ref. [45|] for the QCDF predictions of B u ^ s — > PP, VP and VV decays and 
the comparison with experiment and the theory predictions of pQCD and SCET). This indicates 
that the dynamics of the penguin mechanism in penguin-dominated B decays is well described by 
QCDF. 

As noted in Sec. HI since the charmed quark is not heavy enough and l/m c power corrections 
are so large, a sensible heavy quark expansion in l/m c is no longer applicable to aforementioned 
QCD-inspired approaches. This means that the less sophisticated approaches such as the naive 



factorization or the improved version of factorization such as the generalized factorization [46|, |47 ] 
can serve as a reasonable starting point. However, none of the existing theoretical frameworks 
work well for hadronic D decays as the large l/m c power corrections are nonperturbative in nature 
and cannot be tackled using the factorization approach. Nevertheless, in this work we shall apply 
QCDF to the zeroth order in the power expansion of Aqcd/^c (except for chirally enhanced 
power corrections characterized by the chiral factor r x given in Eq. (|24p which are formally l/m c 
suppressed but numerically very important) to estimate penguin amplitudes in charm decays for 
the following reasons: (i) the expressions of penguin contractions in different approaches are similar 
except that the penguin loop function G(s, x) to be defined below is convoluted with the light-cone 
distribution amplitude of the light meson in QCDF, while it is replaced by G(s,k 2 ) where k 2 is 
the squared momentum carried by the virtual gluon; and (ii) vertex corrections can provide a 
strong phase which is absent in the generalized factorization approach. Hence, we shall work in 



QCDF 42|, |48j] to evaluate perturbative penguin amplitudes, but keep in mind that we employ 



this approach simply for a crude estimate of the penguin contractions. As for power corrections 
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to QCD-penguin exchange amplitude PE and QCD-penguin annihilation amplitude PA, we shall 
adopt a different strategy to deal with long-distance effects due to FSI's, as will be elaborated in 
Sec. UVBl 

To calculate the penguin contributions, we start from the short-distance effective Hamiltonian 



U 



cff 



Gf 
V2 



A p (ciOf + c 2 O p 2 + c 8g 8 g) - XbYl Ci0i 

p=d,s 



i=3 



(21) 



where 



O p 



(pc)v-a( u p)i 



Of 



(p a cp) v _ A (upp a ) 



Oi) V -A1 



03(5) = (mc)v_a^(??) VtA , O m = (u a cp) v _ A J2(qpq a ) 



°8g = -g^m c «cr MI/ (l + 75)G M!/ C 



(22) 



with O^-Oq being the QCD penguin operators and (<jiq2) v±A = 9i7/x(l ± 75)92- The electroweak 
penguin operators are not included in the Hamiltonian as they can be neglected in practice. For 
the Wilson coefficients, we follow 



22] to take ci = 1.22, c 2 



-0.43, c 3 = 0.018, c 4 



-0.046, 



C5 = 0.013, C6 = —0.044 and c 8g = —0.11 evaluated at the scale fx = m c . 
Let us first consider the penguin amplitudes in D — > P\P 2 decays 

Gp 

^| (ifDfpJp 2 ) [bl\ PlP2 
- (ifDfpJp 2 ) K\ Pi p 2 



S PlP2 



pP 



[a 3 (PiP 2 )-a 5 (PiP2)]^ (DPl ' P2) , 
[al{P l P 2 )+r^al{P l P 2 )]X {DP ^ 



(23) 



PA P 



V2 



where C F = (iV c 2 - 1)/(2JV C ) with N c = 3, 

r P (M) 



2m^ 



(24) 
(25) 



x ^ m c (/x)(m 2 +mi)(/t) 

is a chiral factor and 

= (P 2 |jM|o)(P 1 |j;|5) = i/ P2 ( m f,- m | 3 j^( m | 2 ; 

with p c being the center-of-mass momentum of either final state particle. Here we have followed the 
conventional Bauer-Stech-Wirbel definition for the form factor Fq P [l5cj]. The flavor operators a p 
are basically the Wilson coefficients in conjunction with short-distance nonfactorizable corrections 
such as vertex corrections Vi, penguin contractions V% and hard spectator interactions Hi\ 



a 3 (PiP 2 ) 
a p 4 (PiP 2 ) 
a 5 (PiP 2 ) 

<(PlP2) 



C 3 + 

c 4 + 

C5 



Si 

C6 



C5 



c 4 


CfOI s 


N c 




c 3 


CfO-s 


N c 


4tt 


c% 


CfO-s 


N c 


At: 


C5 


C F a s 



[V 3 (P 2 ) + ^-H 3 (P 1 P 2 )}, 
4tt 2 

[V 4 (P 2 ) + —H 4 (P 1 P 2 )]+V P (P 2 ), 
4n 2 



4vr 2 



C6 + ^ ) + ^ -^^(p 2 ) + — i? 6 (P 1 P 2 )] + v p (p 2 



N c 4vr 



iV, 



(26) 
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where the explicit expressions of Vi and Hi can be found in [48j]. The strong phase of the QCD pen- 
guin amplitude arises from vertex corrections and penguin contractions. The order a s corrections 
from penguin contraction read 48J] 



n 



CfOs 
4iTN r 



Cl 



4. m c 2 

o m ~Q G M 2 {Sp) 

6 fi 6 



+ c 3 



8 TTl c 4 

2 m — + 3 ~ Gm 2 (su) - Gm 2 (1] 



P, 



+ (c 4 + Cq) 

C F a s f 

ci 



16 Tfi 

-5- hi Gm 2 { s u) — Gm 2 { s cl) — Gm 2 {s s ) — Gm 2 (1) 



2c 



,eflf 



dx 
1 — x 



M 2 I 



4 , ™c 2 * , 

3 h 7 + r GMj(Sf 



4vriV L v 

+ (C4 + C 6 

where c|g = c§ 9 + C5 , s j = mf / , and 



+ c 3 



8 T m c 4 * - 
gin— + - - Gm 2 {Su) ~ G Al2 {l) 



(27) 



16 771 ^ 

—In— - Gm 2 (s u ) - G M2 {sd) ~ Gm 2 {s s ) - G M2 (1) 



2c 



off 



Gm 2 (s) 



dxG(s, 1 - x)$M a (x), Gm 2 (s) 



dx G(s, 1 — x)$ m2 (x 



with 



G(s,x) 



duu(l — tt)ln[s — u(l — u)x]. 



(28) 



(29) 



Here <3?m 2 (^ > m 2 ) is the twist- 2 (-3) light-cone distribution amplitude for the meson Mi- In the 
generalized factorization approach, the perturbative loop functions Gm 2 (s) and Gm 2 (s) are replaced 
by 

•1 



G(s,kf 



duu(l — u)ln 



k 

u(l - u) — r 

mj. 



(30) 



where k 2 is the virtual gluon's momentum squared. 

The annihilation operators 63 4 in Eq. (|23|) are given by 



N? 



C 3 4 + C5(4 + 4)+^cC64 

c 4 ^i + c 6 Ai 



(31) 



where the annihilation amplitudes A l ^ 2 3 are defined in Ref. 48]. As stressed in [lOj], contributions 
to penguin weak annihilation amplitudes PE and PA from penguin contractions should be included 
in order to ensure renormalization scheme and scale independence. This amounts to the effective 
penguin Wilson coefficients 4g, 1471 ] 



D,eff 
^4,6 



C4,6(m) 



PlM, 



r>,cE 
c 3,5 



C3,5(yw) 



1 



(32) 



2C F 4 ' ovr; ' d '° J ' uv ^ 2C F ' 

Hence, the penguin Wilson coefficients in Eq. (|3ip should be replaced by the effective ones. 

In practical calculations of QCDF, the superscript 'p' can be omitted for 03 and 05. For a 4 6 and 
03^5, their 'p' dependence is dictated by the terms Gm 2 (s p ) and Gm 2 {s p ) defined in Eq. ([28]) . Also 
as explained in [15(, we shall not consider the spectator contributions to Oj because the relevant 
scale Hh « 0.8 GeV is beyond the regime where perturbative QCD is applicable. 
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IV. DIRECT CP VIOLATION 



Tree-level CP violation 



Direct CP asymmetry in hadronic charm decays denned by 

T(D -> f) - T(D -» /) 
T(D^f)+T(D^f) 



a d j P (f) 



(33) 



can occur even at the tree level 5l|]. As emphasized in 15(, the great merit of the topological 
approach is that the magnitude and the relative strong phase of each individual topological tree 
amplitude in charm decays can be extracted from the data. Hence, the estimate of should 
be trustworthy. Larger CP asymmetries can be achieved in those decay modes with interference 

7T°r] and 



between T and C or C and E. For example, a£ is of order (0.7 — 0.8) x 10 3 for D° 



L>+ K + r] (see Table |Tn]). Direct CP violation in D° -> K°K is given by 



(tree) m 
l dir 



2lm(X d X* s ) lm(E* d E s 



1.2 x 10" 



\Ed,E s 



sin 5 ds 



(34) 



I'M \Ed — E s \ \E d 

where 5d s is the strong phase of E s relative to Ed- Substituting the two solutions for Ed and E, 
given in Eq. (|17j) in the above equation yields 



(tree) 
1 dir 



K°K 



-0.7 x 10" 3 
-1.7 x 10" 3 



Solution I , 
Solution II . 



(35) 



. (tree) , ,-.() 

vox comparison, a dir [D 



221. 



K°K ) = 0.69 x 10" 3 is predicted in 
From Table IIIII we see that almost all the predicted tree-level CP asymmetries in 22| are of 
opposite signs to ours. This can be traced back to the phase of the VF-exchange amplitude. For 
CF D -)■ PP decays, its phase is (122 ± 2)° with \ 2 = 0-29 per degree of freedom [Eq. (gj)]. For 
SCS decays, the phases of Ed and E s [see Eq. (fT7|) ] lie in the range of 107° ~ 137°. Therefore, 
the W^-exchange amplitude in this work is always in the second quadrant, while the E amplitude 



m 



22f | lies in the third quadrant because a global fit to all the data of 28 CF and SCS D — > PP 
branching fractions has been performed there with x 2 — 7.3 per degree of freedom. As a result, 



between this work and 



the imaginary part of E in 22| has a sign opposite to ours, and this explains the sign difference 



221 for a 



(tree) 
dir 



B. Penguin-induced CP violation 



Direct CP violation does not occur at the tree level in some of the SCS decays, such as D° — > 
K + K~ and D° — > tt + tt~ . In these two decays, CP asymmetry can only arise from the interference 
between tree and penguin amplitudes denoted by a^ P ^- From Eq. ([7|) we obtain 

4Im[(A d - X s )X* b ] Im[(T* + E* + AP*)(T + E + AP + SP - AP)] n7T 



(t + p), + _N 



|A d — X s 



IT + E + AP 



2 

71 — 



1.2 x 10~ 



P s + PE B + PA S 
T+E+AP 



sin 5-n 



(36) 
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TABLE III: Direct CP asymmetries (in units of 10 3 ) of D — > PP decays, where denotes 
CP asymmetry arising from purely tree amplitudes and 0^°^ from the total amplitude. The first 



(second) entry in D° — > rjrj, rjrf , K + K~ and K°K is for Solution I (II) of Ed and E s [Eq. (fT< 
For QCD-penguin exchange PE, we assume that it is similar to the topological E amplitude [see 
Eq. (I45p ], World averages of experimental measurements are taken from Ref. 0]. For comparison, 



the predicted results of a d * ec ^ and a\^' in [22j are also presented 



(tot) 



Decay Mode 


a£ ee ^(this work) 




a^Vthis work) 

GIT \ 


HIT 1 J 


Expt. 


LJ — 


■> 7T" 1 












0.95 ± 0.04 


0.68 


2.0 ±2.2 


LJ 


*7T° 












0.80 ± 0.04 


0.20 


1 ±48 


LJ 


+ 7T° 






0.82 ±0.03 


-0.33 


0.08 ± 0.04 


-0.55 




u — 


*7T° 






-0.39 ± 0.02 


0.54 


0.01 ±0.02 


1.99 




I J — 


■» r\r\ 




-0.28 ±0.01 


0.28 


-0.58 ±0.02 


0.08 












-0.42 ± 0.02 


0.28 


-0.74 ± 0.02 


0.08 




LJ 


■» r/rj 






0.49 ±0.02 


-0.30 


0.54 ±0.02 


-0.98 












0.38 ±0.02 


-0.30 


0.34 ±0.02 


-0.98 




D°- 


* K + K~ 








-0.42 ±0.01 


-0.50 


-2.3 ± 1.7 
















-0.53 ±0.02 


-0.50 




D°- 




-0.73 


1.11 


-0.63 ±0.01 


1.37 












-1.73 


1.11 


-1.81 ±0.01 


1.37 




D+ - 


-)• 7T" 



















29 ± 29 


D+ - 


->• 7T 


f '/ 




0.35 ±0.06 


-0.54 


-0.74 ± 0.06 


-0.52 


17.4 ± 11.5 a 


D+ - 


-> 7T~ 


f '/ 


1 


-0.21 ±0.04 


0.39 


0.33 ± 0.07 


1.52 


-1.2 ±11.3 a 


D+ - 




-0.07 ±0.06 


-0.14 


-0.39 ± 0.04 


-1.00 


-1.0 ± 5.9 


Df- 


-> IT' 




0.07 ±0.06 


0.14 


0.45 ± 0.03 


1.00 


66 ± 24 


Dt- 




0.01 ±0.11 


0.33 


0.94 ±0.10 


0.72 


266 ± 228 


Df~ 




-0.71 ±0.05 


-0.19 


-0.61 ±0.05 


0.83 


93 ± 152 


Df~ 


-»• if+r/ 


0.36 ±0.04 


-0.41 


-0.28 ±0.12 


-1.78 


60 ± 189 


Data from 


54 


]■ 



where 5 n7T is the strong phase of (P s + PE S + PA s ) nn relative to (T + E + AP) nn . Likewise, 



,(t+p) 

l dir 



-1.2 x 10" 



P d + PE d + PA Q 



T+E - AP 



sin 5, 



(37) 



with S KK being the strong phase of {Pd + PEd + PAd) KK relative to (T + E — AP) KK . Therefore, 
we have the relation 

a d \+ p \K + K-) = -a d t + p) (Tr + Tr-) in SU(3) limit. (38) 
Note that the expression of 0'dir P \ 7r+7T ~) gi yen m fl^ 



,(t+p) 

l dir 



TT IT 



2Im(X d X* s ) Im[(T* + E* + P* d + PE* d + PA* d )(P s + PE S + PA S )] 7 
|A d | 2 

1.2 x 10" 



P s + PE S + PA S 



T + E 



\T + E\l 
sin 5^, 



(39) 
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is a special case of Eq. (|36[) with negligible AP and similarly for a^ p \K + K ). 

Using the input parameters for the light-cone distribution amplitudes of li ght mesons, quark 



masses and decay constants from Refs. 4l|, |52| and form factors from Refs. [3l|, [53|] , we find to the 



leading order expansion in Aqcd/ 7 ™;; in QCDF that 

^) = 0.239 e~ il52 \ (^] = 0.244 e- il54 °, (^\ 



-i35° 



rp I "1 rp j ' I rp 

-*- / 7T7T \ J 7T7T V ^ 

P d\ „„o. -il52° f P s\ _no/in„-il54° ( ^ P \ _ n nnn „-i35 



0.235 e- UM , ^ = 0.240 e" li54 , ^ = 0.009 e"^ 5 . (40) 

Therefore, AP = P^ — P s arising from the difference in the d- and s-loop penguin contractions 
[see Eq. (|26|) ] is very small compared to the tree amplitude. More precisely, it comes from the 
differences between G^^isd) and G Wj k(s s ) and between G nt K(sd) and G^^iss) defined in Eq. (128]) . 
Because of the smallness of AP, we need to rely on SU(3) violation in both T and E amplitudes 
to explain the large disparity in the rates of D° — > K + K~ and ir + ir~ . It is straightforward to find 

Hence, 5 n7T » 5 kk = 170°. From Eqs. (M) and ([37]), we derive af + p) (vr + 7r-) = 6.7 x 10~ 5 
and cl^ p \k + K~) = —4.9 x 10~ 5 . Therefore, QCD-penguin induced CP asymmetries in D° — > 
7r + 7r~, K + K~ are small mainly due to the almost trivial strong phases 5 nn and 5kk- 

For QCD penguin power corrections, we shall consider weak penguin annihilation, namely, 
QCD-penguin exchange PE and QCD-penguin annihilation PA. At the short-distance level, weak 
penguin annihilation contributions are found to be smaller than QCD penguin with the hierarchy 

P > PE > PA. For example, (PE/T) W7r ~ 0.04 and (PA/T) n7T 0.02. In the QCDF approach, 

it can be shown that the short-distance weak annihilation and weak penguin annihilation terms 
are related to each other via (see Eq. (55) of [481 ]) 

^sd = ^2 E SD pA SD = c 4 + c 6 ^SD 
ci c 2 

PE SD = — i£ SD + ^L( ifD f Mi f M2 ) \ C5 Ai + (c 5 + N cC(i )A(\ , (42) 

where A 1 ^ denote annihilation amplitudes induced from (— 2){S— P)®{S+P) 4-quark operator and 
the superscripts V and '/' refer to gluon emission from the initial and final-state quarks, respec- 
tively. For example, the amplitude A3 has the expression A3 oc — 2{M\M2\{uq) s+p ® (qc) s _ p \D) 
with (<?i(/2) g± p = 9i(l ±75)92- Now A3 corresponds to the factorizable contribution 

A{ oc -2(M 1 M 2 \(uq) s+p \0)(0\(qc) s _ p \D) , (43) 

where A\ to the nonfactorizable contribution of A-i. The factorizable term Al was evaluated in 



22| by assuming its dominance by lowest scalar resonances. In [10J, the large- N c counting 

0(N C ) (44) 



(MiM 2 |(to) 5+p ® (uc) s _ P \D) 



{M 1 M 2 \(s a s l 3 - d a d/3) v _ A <g) (uj3C a ) v _ A \D) 

was employed to get the relation PE ~ (2A/ r c Cg ff /c±)E. (The amplitude PE was denoted by P/ 5 i in 
fiol].) However, the major contributions to weak annihilation and weak penguin annihilation from 
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FSI's were not considered in [lol . I22I ]. 4 We would like to stress again that the relations given in 
Eq. (|42p are valid only for short-distance ones. One cannot replace E SU and A SD by the topological 
amplitudes E and A, respectively, extracted from the data. Moreover, annihilation terms at 

SD 

short-distance level are small due to helicity suppression. Typically, QCDF yields E ~ 0.5 and 
j4 SD ~ —0.2 in units of 10 -6 GeV. Comparing with the topological amplitudes extracted in Eq. @, 
it is evident that weak annihilation amplitudes are dominated by long-distance contributions. 

As pointed out in [151 ]. long-distance contributions to SCS decays, for example, D° — > 7r + 7r~, can 
proceed through the weak decay D° — > K + K~ followed by a resonant-like final-state rescattering 
as depicted in Fig. 2 of [3]. It has the same topology as the QCD-penguin exchange topological 
graph PE. Just as the weak annihilation topologies E and A, it is expected that weak penguin 
annihilation will receive sizable long-distance contributions from FSI's as well. Recall that soft 
corrections due to penguin annihilation have been proposed to resolve some problems in hadronic 
B decays, such as the rate deficit problem for penguin-dominated decays and the CP puzzle for 
B° — > K~tt + 48|. Hence, it is plausible to assume that PE is of the same order of magnitude as 
E. For concreteness, we shall take (in units of 10 -6 GeV) 

PE = 1.6e iU5 ° . (45) 

As shown in Table IIIIl we see that the predicted CP violation denoted by a^°^ or a^™^ is at 
most of order 10 -3 in the SM. For Aa^p, the CP asymmetry difference between D° — > K + K~ 
and D° -> tt+it~, we obtain -(0.139 ± 0.004)% and -(0.151 ± 0.004)% for Solutions I and II, 
respectively. It is of interest to notice that the prediction of —0.118% obtained in [22| based on 
a different approach is in agreement with our result. Since in the SM, Aa^ p arises mainly from 
weak penguin annihilation, we can vary the amplitude PE to see how much enhancement we can 
gain. Even with the maximal magnitude \PE\ ~ T and a maximal strong phase relative to T, we 
get Aafjp = —0.27%. This is more than 2a away from the current world average. Hence, if the 
LHCb result for Aa^p is confirmed by further data analysis, it will imply new physics in the charm 
sector. 



C. D+ vr+vr 

It is known that the D + — > 7r + 7r° decay does not receive QCD penguin contributions and its 
direct CP violation vanishes. Nevertheless, it does receive additional contributions from isospin 
violation and electroweak penguin effects. For example, the u-d quark mass difference will induce 
ir-rj-rj' mixing. Hence, the decay D + — > tt + t] (D + — > 7t + t/) followed by the r]-ir (tj'-tt) mixing will 
contribute to the direct CP asymmetry in D + — > tt + tt°. Consider the contribution from 7r°-?78 
mixing 

A(D+ -> 7T+7T ) = -^=\d(T + CV + (7T°|?W|7fe) 2 1 2 A(D + -> 7T+7]) , (46) 



4 Contributions from final-state interactions have the general expression J2 n (MiAl2\Cs\n) (n\Cw\D) , where 
Cs and Cw are strong- and weak-interaction Lagrangians, respectively, and n denotes the physical inter- 
mediate states. 
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where (see e.g., [491 ]) 

f^d ~ m u , 2 2 



4 to, — to ' 



with to = [m u + TO,i)/2. It follows that 

^2 4 m s — m 

There will be also a correction from the rj' meson. Consequently, isospin breaking induced by u 
and d quark mass difference will generate CP asymmetries at tree and loop levels for D + — > 7r + 7r°. 
However, it is suppressed by a factor of (m d — m u )/m s ~ 0.025. Therefore, CP asymmetry in 
D+ — > 7r + 7r induced by isospin breaking and electroweak penguins are negligible. Recently, it was 
argued in 26J that a measurement of nonzero a^p (7r + 7r°) would be a signal for AI = 3/2 new 
physics. 



V. NEW PHYSICS EFFECTS 



Whether the first evidence of CP violation in the charm sector observed by LHCb is consistent 
with the SM or implies NP will require a further analysis of more data and improved theoretical 
understanding. If the improved theoretical estimate of Aa^p remains to be a few per mille and 
the experimental measurement continues to be large with more statistics in the future or if the 
direct CP asymmetry of any of the discussed modes is significantly larger than 10~ 3 , it will be 
clear evidence of physics beyond the SM in the charm sector. Then it will be important to explore 
possible NP scenarios responsible for such large direct CP asymmetries. 



A. Large penguins 

In the wake of the LHCb and CDF measurements of Aa^p, several groups of people have as- 
sumed enhanced hadronic matrix elements through SU(3) breaking [l4| or tZ-spin breaking 21] or 



enhanced penguins via unforeseen QCD or NP effects 16l.ll9l.l24f]. Let us examine the phenomeno- 
logical implications of large penguins irrespective of the origin of enhancement. It turns out that 
the penguin 

i S P = \ Pd + PE d + PA d + P S + PE S + PA S ) « { 2 - 9T C fOT S ° luti0n 1 • (49) 
2 2 V s> \3.2Te 485 for Solution II , V ; 

with maximal strong phase can accommodate the measurement of direct CP asymmetry difference 



between D° — > K + K and 7r + 7r (see Table ITV|) . Indeed, it has been noticed in la, [24| that a 
large penguin of order 3T can explain the observed Aa^p. 5 Some authors argued that penguin 



5 Authors of [3] have introduced an additional phenomenological penguin amplitude Pb in order to ac- 
commodate the measured Aa^p. However, there is no need to make this assumption as the penguin 
amplitudes can be recast to XdPd + A S P S = ^(A^ — X s )(Pd — P s ) — \\b{P<i + Ps)- The last term on the 
right-hand side of the above relation is the so-called Pb in . 
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TABLE IV: Direct CP asymmetries (in units of 10 3 ) of SCS D — > PP decays estimated in 
the scenarios with large penguin contributions and large chromomagnetic dipole operator. The 
parameters £P and c^ p are chosen to fit the data of Aa$ T P : ±£P = 2.9Te i85 ° and c^ p = 0.017e il4 ° 
for Solution I, ^XP = 3.2 Te* 85 and Cgg = 0.012e i14 for Solution II. The number in parentheses 
is for Solution II of E d and E s [Eq. (fTTj) ]. 



Decay Mode Large penguins Large c.d.o. 




matrix elements (or more pertinently "penguin amplitudes") could be substantially enhanced in 
the SM just like the enhancement of the VF-exchange amplitude relative to the naive expectation. 
However, this conjecture needs to be clarified. The topological amplitudes C and E extracted 
from the data are much larger than what expected from naive factorization. This is because they 
receive large l/m c power corrections from FSI's. For example, the topological amplitude E receives 
contributions from the tree amplitude T via final-state rescattering. As shown in Refs. 37J, 



the effects of final-state rescattering via nearby resonances can be quantified. As for the QCD 
penguin, an estimate to the leading-order expansion of Aqcd /fnb i n QCDF is given in Eq. (140 p . 
The l/m c corrections manifested as weak penguin annihilation are even smaller. Nevertheless, weak 
penguin annihilation does receive long-distance contributions from the color-allowed tree amplitude 
followed by final-state rescattering. Since this FSI originates from the tree amplitude, one cannot 
have PE > T. Hence, we believe that a penguin amplitude larger than the tree amplitude in size 
is only possible through the enhancement of new physics. 

Using the above large penguin ^SP as input, the predicted direct CP asymmetries for other 
charm decays are summarized in the second column of Table HVl We see that many modes, such as 
Df — > tt + K°, tt°K + , K + r]' , are expected to yield direct CP asymmetries of a similar magnitude, 
at a few per mille level. 
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B. Large chromomagnetic dipole operator 



Even before the LHCb experiment, the impact of NP on Aajilp had already been investigated 
in 0,[5^]. The unexpected LHCb measurement has inspired many analyses based on a variety of 
NP models. Models for NP effects at tree level include flavor-changing coupling of a SM Z boson 
3)0], flavor-changing neutral currents induced by a leptophobic massive Z' boson [lllQ], two 
Higgs-doublet model [19[], color-singlet scalar model [131 ]. color-sextet scalar (i.e., diquark scalar) 



model 



0,E] 



color-octet scalar model 



19f ] and fourth generation model 12], [2]]]. Models with 



NP in QCD penguins at the loop level have been constructed as well, inclu ding new fermion and 
scalar fields [19f] and the chirally enhanced chromomagnetic dipole operator [la ]. 

The NP models are highly constrained by D°-D mixing, mixing and CP violation in 

the kaon system characterized by the parameter e'/e 0]. Many of the tree-level NP models are 
either ruled out or in tension with experiments [19j] . As stressed in [18|] , a large NP contribution to 
the AC = 1 chromomagnetic dipole operator is the best candidate to explain the LHCb and CDF 
results as it is least constrained by all current data in flavor physics. Although the chromomagnetic 
dipole operator Osg is suppressed by the charm Yukawa coupling, the hadronic matrix element reads 



(M 1 M 2 |0 8s |L») = ^(M 1 M 2 \u a ^(l + 75)^^ q l7 ^ qs \D) 



(50) 



where k 2 is the square of momentum transfer of the gluon and is of order m 2 . After applying the 
equation of motion, we see that the matrix element {MiM2\Os g \D) is independent of m c ; that is, 
it is enhanced by a factor of v/m c , where v is the vacuum expectation value of the Higgs field 
arising from the structure of the gauge-invariant dimension-six operator 18j. On the other hand, 
the D°-D° mixing induced by 0$ g is suppressed by a factor of m 2 /v 2 . Of course, we need NP to 
enhance the Wilson coefficient c$ g and to induce a sizable imaginary part. This can be realized in 
the supersymmetric models where the gluino-squark loop contributes a major part of c% g floT ], the 
disoriented A terms and split families are the sources of flavor violation 18|] , or the flavor structure 
of the trilinear scalar couplings is related to the structure of the Yukawa couplings via approximate 
flavor symmetries 251 ]. 

To demonstrate the NP effects, we consider the NP penguin amplitude P NP induced by the 
chromomagnetic operator which has the expression 



3 NP 
P1P2 



^[aT(PiP2)+r?ar(PiP2)]X( DP ^ 



(51) 



with 



«4 NP (^2) 



NP 
a 6 



rf p (P2) 



(PP 2 ) = p 6 np (P 



-2c 



np 



-2c 



NP 



CfOs 
4itN c 

CfOCs 



d.r 



1 



■$P 2 (x) 



(52) 



AttN c 

See Eqs. (|26l) and (|27|) for a derivation. Then we add the NP amplitude P NP to the penguin 
amplitudes \ p P p in Table [I] where summation over p = d, s is understood; that is, the penguin 
amplitudes X p P p are replaced by X p P p + P NP . As an example of illustration, we shall take c$ g p « 
0.012e ?14 ° which fits to the data of Aa^p. The calculated CP asymmetries for the other modes are 
listed in the last column of Table HVl Since the decay D° — > K°K° does not receive QCD penguin 
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contribution, it is not affected by the chromomagnetic dipole operator in NP. It is interesting to 
notice that while a large chromomagnetic dipole operator leads to a large direct CP asymmetry 
for D° -> 7r 7r ,vr 7?, the predicted CP violation for D° -> vr ?/, D + -> k + t/,K + K° and D+ -> 
ir + K°,K + r]' is much smaller than that in the large penguin scenario. Therefore, measurements of 
the CP asymmetries of the above-mentioned modes will enable us to discriminate between the two 
different NP scenarios. 



VI. DISCUSSIONS AND CONCLUSIONS 

In this work we have examined various sources responsible for the seemingly large SU(3) breaking 
effects in the D° — > K + K~ and D° — > 7r + 7r~ decays. We considered three cases: (i) SU(3) 
symmetry holds for T and E amplitudes. Then a sizable AP (the difference of s- and <i-quark 
penguin contractions) of the same order of magnitude as T is needed to explain the data. This 
is also true if the symmetry breaking of T and E amplitudes follows the pattern given by by Eq. 
(fTTj) . (ii) SU(3) symmetry holds for E amplitudes, but SU(3) violation due to decay constants, 
meson masses and form factors is taken into account in T amplitudes so that T KK /T V ^ ~ 1.32 
as inferred from the factorization approach. This leads to \AP/T\ ~ 0.15. (iii) The large rate 
difference between D° — > K + K~ and D° — > tt + it~ is entirely accounted for by SU(3) violation in 
T and E amplitudes and hence AP can be neglected. Owing to the observation of D° — > K°K 
through ^-exchange and penguin annihilation diagrams and the smallness of AP theoretically, 
we have argued in this work that the last scenario is preferred and fixed the SU(3) breaking in 
the W-exchange amplitudes from the following four modes: K + K~, 7r + 7r~, 7r°7r° and K°K . Our 
results are summarized in Table HU where we have shown branching fractions of SCS D — > PP 
decays and elaborated on SU(3) breaking effects for each mode. 

Since the magnitude and the phase of topological color-suppressed tree amplitude and weak 
annihilation amplitudes which arise mainly from final-state rescattering can be extracted from 
the data in the diagrammatic approach, direct CP asymmetry aj^ at tree level can be reliably 
estimated. We predict that a { ^ e) (D° -)• K°K°) ranges from -0.73 x 10" 3 to -1.73 x 10~ 3 , 



depending on the solution for E amplitudes. A recent similar study in 22j found opposite signs of 



a ^dir^ t° ours for most of the SCS D — > PP decays. This is ascribed to the phase of E amplitudes: 
It is in the second quadrant in our work while in the third quadrant in 22|. As for the decay 
D + — > ir + ir , it does receive corrections from isospin violation due to the u-d quark mass difference 
and electroweak penguins. However, the induced CP asymmetry is negligible because, for example, 
the isospin- violating effect is suppressed by a factor of [ma — m u )/m s . 

Using QCD factorization as a guideline, the direct CP asymmetries of both D° — > K + K~ and 
D° — > 7T + 7r _ are at a few xl0~ 4 level. This is seen to be largely due to the trivial relative strong 
phase between the QCD penguin amplitude and the tree-level amplitudes. For QCD penguin power 
corrections, the short-distance contributions to weak penguin annihilation diagrams PE and PA 
are small, but PE receives long-distance final-state contributions from rescattering through nearby 
resonances which have the same topology as the VF-exchange diagram. It is thus natural to assume 
that PE is of the same order of magnitude as E. We conclude that the CP asymmetry difference 
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Aa^p between D° -> K+K' and D° -> vr+vr" is about -(0.139 ± 0.004)% and -(0.151 ± 0.004}% 
for the two solutions of E amplitudes, respectively. A similar prediction was also obtained in 22f | . 

If Aa^'p continues to be large with more statistics in the future, it will be clear evidence of 
physics beyond the standard model in the charm sector. Considering two possibilities of new physics 
effects, namely, large penguins and large chromomagnetic dipole operator, we have studied their 
phenomenological implications in the SCS charmed meson decays. We point out that the CP asym- 
metries several modes, such as D — > —7- 7T + 7/, K+K and D+ -> tt+K , K+rf, 
allow us to discriminate between different new physics scenarios. 
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